Nonlocality and fluctuations near the optical analog of a sonic horizon 
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We consider the behavior of fluctuations near the sonic horizon and the role of the nonlocality 
of interaction (nonlinearity) on their regularization. The nonlocality dominates if its characteristic 
length scale is larger than the regularization length. The influence of nonlocality may be important 
in the current experiments on the transonic flow in Kerr nonlinear media. Experimental conditions, 
under which the observation of straddled fluctuations can be observed, are discussed. 
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I. INTRODUCTION. 

A possibility of creating a sonic horizon in a 
transonically accelerating flow^, mimicking either 
black or white hole horizons, inspired proposals 
and attempts to realize them. Quite a number 
of theoretical models realizing such artificial black 
(white) holes has been put forward 2-7 . Exper- 
imentally a white-hole horizon was observed in 
optical fibers^ where the probe light was back- 
reflected from a moving soliton, and a black-hole 
horizon was observed in a Bosc-Einstcin conden- 
sate (BEC) system^. A radiation in optical fibers 
was observed^, which is a promising contender 
to being analogous to Hawking radiation in a 
table-top experiment. A "horizon physics" is cur- 
rently intensively studied also in the surface water 
wavesii~— . The whole field of analog gravity has 
been recently reviewed in Ref. [bH 

Propagation of the coherent laser beam in Kerr 
nonlinear media is in its many aspects analogous to 
a flow of fluid. Properties of such "luminous fluid" , 
in particular formation of the so called " dispersive 
shock waves" were studied theoretically-^— and 
experimentallyi^r— . Similar approach appeared to 
be very successful in studying temporal dynam- 
ics of tunneling 2 ^—. This analogy prompted a 
proposal 2 ^ to create an optical analog of the Laval 
nozzle, in which the luminous fluid is accelerated to 
a supersonic velocity and a sonic horizon is created. 
Our recent experiment 2 ^, realizing the proposal 2 ^, 
demonstrated a possibility of transonic accelera- 
tion of a luminous fluid in an all-optical analog 
of the Laval nozzle. A laser beam is propagat- 
ing in a specially profiled cavity filled with ethanol 
with an addition of iodine. The beam disperses 
and also deviates from the straight line (along the 
z axis) propagation so that the beam angle with 
respect to the axis mimics the flow velocity and 
its variation mimics acceleration. The defocusing 
nonlinearity in this experiment appears due to a lo- 
cal temperature variation under the action of the 



laser beam. Such mechanism directly leads to non- 
locality of the Kerr nonlinearity (see review 2 ^) . It 
was demonstrated both theoretically and experi- 
mentally in Ref. that the nonlocality can be 
described by a nonlocal response function with a 
finite radius. The nonlocality in Kerr nonlinear 
media allows one to control stability and mobility 
of solitons (bright or dark, depending on the sign 
of the nonlinearity) as shown theoretically 2 ^— and 
confirmed experimentally. 31 In Bose-Einstein con- 
densates the nonlocality may become of impor- 
tance at high enough densities^. It is also typi- 
cal of the dipolar bosonic quantum gases^ 3 - and of 
nematic liquids 3 ^. 

However the issue of nonlocal nonlinearity has 
not yet been addressed in the context of the hori- 



zon physics. A number of papers 



5.6.35-40 



present 



detailed analysis of fluctuations and Hawking radi- 
ation in transonically accelerated fluid accounting 
for local nonlinearity only, which is a quite rea- 
sonable approximation in BEC at least for very 
short-range interactions. The same approximation 
in Kerr nonlinear media may be less harmless and 
an analysis of the role of nonlocal nonlinearity in 
dynamics of fluctuations and formation of Hawk- 
ing radiation near the Mach horizon is necessary. 
The nonlocality introduces a new length a (see Ref. 
I27I ) so that we may expect that fluctuations with 
shorter lengths will be smeared out. The interplay 
of the nonlocality length a with other lengths char- 
acterizing the problem will be discussed in what 
follows and conditions for observation of analog 
Hawking radiation in nonlocal Kerr nonlinear me- 
dia will be formulated. 



II. NONLOCAL NONLINEARITY. 

We consider here the Nonlinear Schrodinger 
(NLS) equation in 1 + 1 dimensions, 



id z ^(x, z) 
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where ^ is a paraxial amplitude of light propagat- 
ing along the z axis, which plays the part of time. 
(3 is the wave vector of the laser beam, which plays 
now the role of mass. R is a linear integral nonlocal 
operator, 



dx'dz'R{x - x, z' - z)\^{x',z')\ 2 (2) 



The kernel of this operator is normalized to one 
and is assumed to be not singular and character- 
ized by a finite length scale a. For example it may 
be chosen in the form 
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In Kerr nonlinear optics the nonlocality depends 
obviously both on x (" coordinate" ) and z (" time" ) . 
The dependence is expected to be symmetric. The 
Fourier transform of §5§ reads 



R{k x , k z 
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(1 +<T 2 k 2 x + (T 2 fc2)3/2- 



It tends to one in the long wave limit k x ,k z — > 
0, and to zero in the short wave limit k x ,k z — > 
oo. These generic features (see, e.g. discussion in 
Ref. Elh will be important for the analysis to be 
presented in the following. 

The Madelung transformation — fe~ lip allows 
one to represent the NLS equation for a complex 
function as two hydrodynamic equations 



d z p + d x (pv) = 0, 



\d x v 2 = -id x 



U, 



qu 
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(5) 



for two real functions. Here 



2/3 / ' 



These two equations describe an equivalent lumi- 
nous fluid, where the light intensity p = f 2 plays 
the role of the density and f3v = —d x tp defines the 
flow velocity v in the x direction. 

Our aim here will be to analyze the behavior of 
fluctuations near the Mach horizon of the transoni- 
cally accelerating luminous fluid. This analysis will 
be quite analogous to the one described in detail in 
Ref. [38| (references to some earlier papers can also 
be found there). The only difference is due to the 
nonlocality term in the Euler equation ([5]). That 
is why only principal steps will be outlined below, 
which are necessary to introduce the notations and 
arrive at the result. 



We consider fluctuations on the background of 
a given stationary flow profile described by a func- 
tion "to = foe~ llpo . They can appear spon- 
taneously or induced artificially by experimental 
means. 24 Linearized equations for these fluctua- 
tions are deduced from 0} and ©, 

(d z + v d x )x - -^—d x {p Q d x i) = 0, 
PPo 

(d z + v d x )€, + ij^d x (pod x x) ~ g(R-Pox) = 0. 

where the density 5p(x, z) and velocity Sv(x, z) 
fluctuations arc defined by the relations 



Sv(x,z) = -~d x £(x,z), 
Sp(x,z) = p (x)x(x,z), 



(7) 



and p = f$, v = -(1/P)d x ip . 

Equations resulting from linearization around 
a known solution are usually called modulation 
equations, which may or may not lead to a modu- 
lation instability. Such an instability appears typ- 
ically in the case of negative g < 0, i.e. focusing 
nonlinearity, whereas we deal here with positive, 
g > 0, i.e. defocusing nonlinearity. Refs. I4lll42l 
present a general discussion in the context of the 
Kerr nonlinear media. In the simple case of a con- 
stant background density po and velocity vq we 
readily obtain the spectrum of the fluctuations in 
the form 



(k z - v k x f 



gpo 
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k 2 l 2 



R(k Xl k z ) 



, (8) 



which corresponds to the Bogoliubov excitation 
spectrum in the case of nonlocal interaction. Here 
l 2 t = 1/(2/3(7 Po) is the nonlinearity length (healing 
length in BEC). If the nonlocality kernel is chosen 
in the form ([3]) then the long wave limit in (|8]) holds 
at k x a 1 and k z a 1 and the sound velocity 
keeps its standard form s 2 = gpo/P- We assume 
the nontrivial situation when a > l n . Then the 
quartic dependence 

k 4 

(k z - v k x ) 2 = ^ 

becomes dominant under the condition that k x z > 
(lla 3 )- 1 / 5 , which may be fulfilled even at l n k x _ z < 
1. 



III. REGULARIZATION DUE TO THE 
NONLOCALITY. 

In order to consider the fluctuations close to the 
Mach horizon we will use the relations 



v = s(l + ax), 



gpo(x) 



ax) 
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describing a spatially accelerating flow. Here s is 
the sound velocity at x = 0, i.e. at the throat of 



where the terms 0(a 2 x 2 ) are omitted. These equa- 
tions can be solved in the Fourier space to within 
the terms 0(a/k x ) (see Ref. [38h - The Fourier 



the Laval nozzle. Then Eqs. ([6]) take the form 



(9) 



transformed first equation Q is solved with re- 
spect to which is then substituted into the sec- 
ond equation. As a result we get 



d zX + 3(1 + ax)d xX + - R [ad x - dl] £ = 0, 



— [-ad x + dl] X ~ /3s 2 Rx + ps 2 aR(x X ) + [d, + 5(1 + ax)d x ]t = 0. 



dk x In Xk 



-^(2iak x + k x ) + (iak x + k x )R(k x , v ) — (y — k x ) — iak x 
ak x {2v - [2 + R(k x , v)]k x - iaR(k x , v)} 

I 



(10) 



where v — k z /s. 

The integration of the r.h.s. of Eq. (fTU|) . al- 
though possible, may result in very cumbersome 
expression. That is why we consider here two lim- 
its. If ak X}Z <C 1, then R{k Xl v) ~ 1 and we get 

dk m hxxk{k x ,v) = 



{ ~ ak x (2„-3k x - la) • 

It means that we return to the situation of the 
local nonlinearity. As shown in Ref. 38 there is a 
singular real space solution of Eq. (|9|) 
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where 



Xs{x,v) oc x 
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Finally it results in the ^-dependent Hawking tem- 
perature 



T H (u)=T H (0)/(l+-l 2 y) 



where 



Th(0) 



3hs 
Airks 



This solution holds for the distances from the Mach 
horizon satisfying the condition min{l/V, 1 /a} S> 
|a;|. In the absence of nonlocality the other con- 
dition is |af| 3> l r where l r — Zn/Gn Q 1/ ' 3 is the 
regularization length^ (see also^) . As we will see 
below the nonlocality also leads to a regulariza- 
tion which means that the final condition reads 
| a; | max{cr, l r }. 

We also have to consider the limit k x z a ^> 1. It 
produces nontrivial results only if a > l n . Other- 
wise the problem can be reduced to the local one. 
In the limit k x a 3> 1 when R(k x , i/) < 1 equation 
([TO]) becomes 

dk m lnxk(k x ,v) = 



" (2iakl + k x ) - v 2 + 2vk x - k' x - iak x 



■ 2 



2ak x (v — kx) 



(14) 



Carrying out the procedure, as outlined in Ref. l38l 
the real space density fluctuations for a given v are 
described by the function 

where 



7i 



2a' 



12. 2 2 4a 
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The behavior of the density fluctuations can be 
readily deduced from (fT5j) . At v <C l~ x and a 3> 

|x| > l n we get 

^k+\kllv-^llv 2 k+\llk 2 -^kHlv--l-llk\ 
2a 2 " 4a " 4 " 8a ™ 12a " 



and 



A = 



-ikz 



dkk^ (k 



-ik z z 



2a 



dyy 11+ ~< 2 e iv 



(16) 



where integration variable y = [x + l/(2a)]fc is 
used. It means that the function \ tends to a con- 
stant (since a\x\ <C 1) and the singular solution 
(fT2|) becomes regular at small enough x. This reg- 
ularization is not violated at even smaller distances 
\x\ <C l n and may be obtained in the way similar 
to that used in Ref. [H. 



IV. DISCUSSION AND SUMMARY. 

As discussed in Introduction, the non-locality of 
the nonlinearity may play an important part in ex- 
periments designed to observe Hawking-type clas- 
sical and quantum fluctuations near the Mach hori- 
zon formed by a transonic flow of luminous fluid. 
The above derivation shows that the non-locality 
should be taken into account if its characteris- 
tic length scale exceeds the regularization length, 
a > l r . The largest of these determines the point , 
below which the singular fluctuation eigenfunction 
([12} is regularized. Consequently, if a > l r it is 
the length scale a, characterizing the nonlocality, 
which sets the lower bound for the distance \x\ such 
that Hawking radiation or straddled fluctuations 
appear. 

At the same time it should be taken into ac- 
count that the derivation of fll~2|) holds for distances 
from the horizon that are not too large, and in par- 
ticular shorter than 1/a (the inverse acceleration 
rate) and 1/v (the length scale of a fluctuation in 
the z direction). The latter sets an upper limit 
for the distance from the horizon where a fluctu- 
ation may develop (see e.g. Ref. 24 for a dis- 
cussion of a Kerr nonlinear medium). A window 
min{l/i/, 1/a} > \x\ > max{/ r ,a} in real space 
must therefore exist in order for fluctuations to be 
observed experimentally. In a typical experiment, 
such as the one described ir*2£, the nonlinear length 
is on the order of several tens of microns, while 
the non-locality length is on the order of 1 mm21. 
This means that 1/v and 1/a must both be at 
least on the order of a few millimeters. 1/v can 



be as large as the size of the experimental appa- 
ratus (a few centimeters), and is not expected to 
be the limiting constraint. On the other hand, the 
experimental setup must be designed in such a way 
that the acceleration of the luminous fluid should 
occur over a distance of several millimeters in or- 
der that the required window for the fluctuations 
exists. Changing the stationary profile allows one 
to change the spatial acceleration a and therefore 
l r . This offers the possibility by increasing a to 
turn to the regime where a is the length scale sep- 
arating the regular from the singular behavior for 
all a > a c = (l n /a) 2 /a. In case when the quan- 
tum potential U qu can be neglected, i.e. we can 
assume that l n = 0, it is the nonlocality regulariz- 
ing the fluctuations for |x| < a. The relationships 
between these different length scales are schemat- 
ically shown in Fig. [TJ 




FIG. 1: (Color online) A sketch of an all-optical Laval 
nozzle demonstrating the relative length scales which 
determine the window in real space where straddled 
fluctuations may be generated in an experiment, a is 
the non-locality length (typically 1 mm); l r is the reg- 
ularization length, determined by the nonlinear length 
and the acceleration rate (typically a few tenths of a 
mm); a is the acceleration rate (typically > 1mm" 1 ); 
v is the characteristic "frequency" of the fluctuations, 
and should be such that 1/v be of the order of the size 
of the experimental test cell (a few centimeters). 

To summarize, we considered the effect of the 
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nonlocality of the nonlinear response on the feasi- 
bility of observing, in an actual experiment, fluc- 
tuations near a " sonic" horizon of an accelerating 
luminous fluid that are analogs of Hawking radi- 
ation. We find that the non-locality sets a lower 
bound on the distance from the horizon where the 
fluctuations may develop, and this, in turn, sets an 
upper bound on the acceleration rate of the fluid. 



With proper design of the experimental apparatus, 
there appears to be a large enough window for ob- 
serving such fluctuations. 
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